Floating-Point-Roundoff Accumulation 
in Digital-Filter Realizations 

By I. W. SANDBERG 

(Manuscript received June 20, 1967) 

In this paper, several results are presented concerning the effects of 
roundoff in the floating-point realization of a general discrete filter governed 
ideally by a stable difference equation oj the form 



IV 



n = Z) hXn-h ~ 2 «*«>»-* n ^ N (1) 



fe-i 



in which [w„] and {x„\ are output and input sequences, respectively. 

In particular, for a large class of filters it is proved that there is a func- 
tion f(K) with f{K) — * as K — » °o and a constant c, both dependent 
on the b k , the a k , the order in which the products on the right side of (1) 
are summed in the machine, and t, the number of bits allotted to the mantissa, 
such that 

(e) K ^ c(y) K + f(K) 

for all K ^ N, in which, with [y n ] the computed output sequence of the 
realized filter, 



and 






Bounds on f(K) and c are given that are not difficult to evaluate, and which, 
in many realistic cases, are informative. For example, for the second-order 
bandpass filter: 

w n = x n — aiiv n -! — a- 2 w n - 2 , n ^ 2 (2) 

with a, and a 2 chosen so that its poles are at approximately ± 4-5° and 
at distance approximately (but not less than) 0.001 from the unit circle, 

1775 
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we find that c, an upper bound on the "asymptotic output error-to-signal 
ratio", is not greater than 0.58 X 10'*, assuming that t = 27, that the 
terms on the right side of (2) are summed in the machine in the order in- 
dicated (from right to left), and that the x n in (2) are machine numbers. 
If the x n are not machine numbers, and hence must be quantized before 
processing, then c ^ 0.76 X 10~ A . 

In addition to error bounds, an inequality is derived which, if satisfied, 
rules out certain types of generally undesirable behavior such as self- 
sustained output limit cycles due to roundoff effects. This inequality 
is satisfied for the example described above. 

I. INTRODUCTION 

The difference equation 

M N 

w n = 2 &frC B _i - 2 a k w n - k , n ^ N (1) 

with M ^ N defines the behavior of a general time-invariant discrete 
filter which acts on an input sequence x , x x , x 2 , ■ ■ • to produce an 
output sequence w N , w lV+l , w lV+2 , ■ • • that depends on the starting 
values Wo , Wi , • • • , %_i . 

There is a vast literature concerned with techniques for designing 
discrete filters [i.e., for determining the a k and the b k in (1)] to meet 
specifications of various types (see, for example, Refs. 1, 2, and 3), and 
a good deal of material is available on the subject of roundoff effects 
in fixed-point realizations of discrete filters (see, for instance, Refs. 4 
and 5). In this paper, we derive some bounds on a meaningful measure 
of the overall effect of roundoff errors for discrete filters realized as 
digital filters on a machine employing floating-point arithmetic oper- 
ations. This type of realization, as opposed to the fixed-point kind, is 
of particular importance in connection with, for example, digital com- 
puter simulations of systems, as a result of the large dynamic range 
afforded by the floating-point mode. 

There are basic differences concerning fixed-point and floating-point 
error estimation problems which stem from the fact that the modulus 
of every individual arithmetic error in the fixed-point mode is bounded 
by a constant determined by the machine, whereas the maximum 
modulus of the error in forming, for example, the floating-point sum 
of two floating-point numbers is proportional to the magnitude of the 
true sum. For this reason, the approach* presented here, as well as the 



♦The approach can be extended in several different directions. For example, 
it can be used to obtain statistical error estimates based on the assumption that 
each roundoff error is an independent random variable. 
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character of the results, are quite different from those of earlier 
writers concerned with fixed-point realizations. 

In addition to error bounds, an inequality is derived which, if 
satisfied, rules out certain types of generally undesirable behavior 
such as self-sustained output limit cycles due to roundoff effects. 

II. ASSUMPTIONS AND RESULTS 

2.1 Assumptions 

It is assumed that : 

(i) each machine number q is equal to sgn (q) a 2" In which 
the exponent b is an integer, and a, the mantissa, is a t-b'it number 
contained in [i, 1] or [f, 1] U {0|; 

(ii) the range of values of b is adequate to ensure that all computed 
numbers lie within the permissible range; 

(in) the machine operations of addition and multiplication are per- 
formed in accordance with standard rounding conventions* (described, 
for example, by Wilkinson ) ; and 

(iv) the coefficients a k and b k in (1) are machine numbers. t 

2.2 Results: x n Machine Numbers 

It is assumed throughout Section 2.2 that the x„ of (1) are floating- 
point machine numbers. 

If the discrete filter (1) is realized on a floating-point machine, then 

y„ = Slit &**-* - t a k y n -) , n^N (2) 

\A=0 ft-1 ' 

in which the y n arc approximations to the infinite precision numbers 
w n , and //(2 — S) denotes the machine number corresponding to 
(S — S) with the understanding that the floating-point numbers cor- 
responding to the products &*s»-* and a k y n . k are to be machine-added 
in some specified order. 
Let 

D(z) ^ 1 + Z a k z-\ (3) 



♦That is, conventions for which the first two equations of Section III am 

satisfied. . . , „ 

fit is certainly true that preliminary design considerations may lead to coeili- 
cients that arc not machine numbers, and one may then be interested also in 
the overall effect of approximating Hie coefficients by machine numbers. That 
problem also can be treated with the approach used here. 
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let 



<«>* = GfTT 1 1 & I 2 )' 



for every sequence {q k } and all K ^ 0, and let e„ denote the nth error 
(y» - u> n ) for n ^ 0. 

Our first result (all proofs are given in Section III) is as follows. 
If D(z) ?* for | z | ^ 1 [i.e., if the discrete filter (1) is stable], then 

(e) A - ^ max | D(e iu )~ x \ [ K £ I V* \j 

+ 2-(E I h | ft) max | D(e'T l I GrTT ^ I *" 

+ 2-'(Z I a k | a,) max [ D(e ia y l I <Z/>* (4) 

for all K ^ N, in which, with y n = w n = for n < 0, 



v 



. = E a*(2/»-* - <»n-*) n = 0,1,2, ■-■ ,(N -1) 



*-0 



and the a* and & are easily evaluated nonnegative numbers which 
depend on the order in which the products in (2) are summed. 

Since the first term on the right side of (4), which arises as a result 
of the possibility of differences in the starting values, approaches zero 
as K — > oo, we see that, after a reasonable number of evaluations of 
the successive y n , (e) K is bounded essentially by a constant times the 
root-mean-squared value of the input sequence, plus another constant 
times the root-mean-squared value of the output sequence. 

In order to determine the a k and /3* , we draw a signal-flow graph 
that indicates the ordering of the operations that would be used to 
compute 



/ u N 

\jfe-0 *=1 



(5) 



if x n and y n were unity for all n. This graph is to contain an input 
node with input b' k for each b k ^ 0, an input node with input a' k for 
each a k 9± 0, no other input nodes, and a single output node 6 which 
is associated with 

£ b' k - Z a' k . 
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All other nodes represent an addition or subtraction of two signals 
to produce a third signal. Exactly one branch is connected to each 
of the input nodes and to the output node. We assign the value p to 
all of the branch transmissions with the exception of those branches, 
if any, which terminate on an input b' k or a k for which b k or a k , respec- 
tively, is equal to unity. These branches are assigned unity transmission. 
Then, by inspection, we evaluate the signal at 6, which must clearly 
be of the form 

Zwp" (4) + E«Jp m,m (6) 

in which ippQc) and <p a (k) are positive-integer valued functions. In 
terms of these functions* 

ft = (1.06)^(fc) 

a t = (im)<p a (k). 

For example, if the right side of (2) is computed as the floating-point 
difference of the machine sums 

fl(b x n + &,*„_, + • • • + b„x n - M ) 

and 

flfaiyn-i + a,y n - 2 + • • ■ + a w y„- N ), 

each obtained by performing machine summations in the order indicated 
(from left to right), if all of the b k and a k are nonzero and not unity, 
and if M ^ 1 and N ^ 2, then the relevant flow graph is shown in 
Fig. 1, from which it follows that 

A, = (1.06)(M + 2) 

ft = (1.06)(M + 2) 

ft = (1.06) (3 + M - k); k = 2, 3, • • • , M 

ai = (1.06)(iV+ 1) 

a 2 = (1.06) (AT + 1) 

a k = (1.06)(3 + N - k); k = 3, 4, • • • , N. 

The bound (4), although revealing, requires a lcnowledge of both 
(x) K and (y) K and is, therefore, not as explicit as we would like. 

* We are assuming here only that max \q<g(k)\2'' < 0.1 and max \<p a (k)\2~' 
< 0.1. Also if tpfi(k) = 1, then we can take /3 t = 1, and similarly for <p„(fc). 
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Fig. 1 — Flow graph for the example. 

For the important case in which b ^ and N(z) = z2lL b k z~ k 9^ 
for I z I ^ 1 (i.e., for the minimum-phase filter case) we prove that 
if the filter (1) is stable and if 



min I N(e iu ) I > 2-' £ I 6* I ft , 

0§wS2tt 4 = 



(7) 



then there exists a constant c, independent of K, and a function j(K) 
with the property that }(K) — > as K — > 00 such that 



(e) K ^ c(y) K + KK) 
for all K ^ N. Moreover, it is proved that 



(8) 



^ 2 - ' max I D(e ia )' 



£ I «* I «* + E I M & 



max I D(e ia )/N(e iu ) | + max | i\T(e m '") _1 [ 2 * £ I «* I «* 



1 - 2"' £ I b t I ft max | Ntf")'* I 



(9) 



and 



/(*) =g max I ZW I (-^rj g I * I 2 ) + m^ I WT' I 2 
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(10) 



, u v max j N(e"°) ' | (t^-XT S I 9« T ) 

• (Z I ^ I ft) s \±±±^ L 

1-2- 1 E I 6. I A max I W")" I 

*• = II u 

for all K ^ A r , in which, with a = 1 and x n = y„ = for n < 0.. 

.V M 

*=0 4=0 

for n = 0, 1, 2, • • ■ , (A' - 1). 

Since <//) A - is the root-mean-squared value of the computed output, 
and since /(/v) -» fairly rapidly as K — » «>, we may interpret the 
smallest value of c for which (8) is satisfied (for all input sequences) 
as an "output error-to-signal ratio" of the realized digital filter. Note 
that the bound (9) on c is not difficult to evaluate. 

2.2.1. Stability in the Presence of Roundoff 

If roundoff effects are ignored, it is well known that the discrete 
filter is stable in several different senses of the word if D(z) ^ for 
\z\ ^ 1. In Section III it is proved that, with roundoff effects taken 
into account, the digital filter is stable in the sense that there is a 
constant c, and a function /,(/v), with j x (K) independent of the values 
of .r„ for n :> N and f^K) — > as K -» co , such that 

(ll) K ^ c,(x) K + UK) (11) 

for all A' ^ N, provided that D (z) 3* for | 2 | ^ 1, and 



mm 



W) I > 2"' Z I <** I «i • ( 12 > 



Roughly speaking, inequality (12) is satisfied if the damping of the 
infinite precision counterpart of the digital filter is sufficiently large 
relative to the number of bits allotted to the mantissa. Stability in 
the sense of (11) rules out, for example, the possibility, due to roundoff 
effects, of a limit-cycle response to a zero input sequence or to an 
input sequence [x n ] that approaches zero as n — > <*> .* 

* There are simple examples which illustrate that instability may result with 
D(z) 9± for \z I ^ 1 if ( 12) is not satisfied. For instance, suppose that each machine 
number is represented in the form ( — »m.2° + ?/n2 _1 + m&~ % + ■ ■ ■ + m,2 - ')2'' 
with the nij zeros or ones, and t > 1. Let 

«•„ = (1 - 2-')»' n _, + (1 - 2-')2-'w„_2 for n ^ 2, with 10 = "'i - — 1. 

Then//[(1 - 2-')u>,] = -(1 - 2-'), ffl(l - 2-')2-'«'o] = -(1 - 2-')2"', and 
fl[-(l - 2~') - (1 - 2"')2 _ '] = — 1, which shows that the computed approximation 
y„ to «■„ satisfies w„ = — 1 for all n ^ 0. This example is a slight modification of one 
suggested by S. Darlington. 
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2.3 A Result Concerning the Overall Effect of Input Quantization Errors 

In many applications the sequence {x n } of (1) is obtained by quantiz- 
ing an input sequence {x n } [i.e., by replacing each x n with the machine 
number (or one of the possibly two machine numbers) of closest value]. 
The infinite precision response w N , w w+1 , •••to the sequence {x N \ 
satisfies 

M N 

u\ = X b k x n - k — X] a>cW„-k , n ^ N (13) 

*=0 4=1 

with w , Wi , • • • , «)#_, some set of starting values. Let w N , w>. v+1 , • ■ • 
be defined by (1) with w n = w n for n = 0, 1, 2, • • • , (N — 1). It is 
clear that (y — w) K , the root-mean-squared value of the difference 
of the computed output and the infinite precision response to {x„}, 
satisfies 

(y - w)k ^ (y - w)k + (w - w) K . (14) 

Bounds on the first term on the right side of (14) are given in Sec- 
tion 2.2. In Section III it is proved that if both N(z) and D(z) have 
no zeros on or outside the unit circle, b ?* 0, and 

min | N(e ia ) \ > 2~' £ | b k \ k , 

u ft = 

then* there is a constant c 2 and a function f 2 (K) such that j 2 (K) — > 
as K — > oo , and 

(w - w) K =g cMk + UK) (15) 

for all A" ^ N. It is proved also that 

M 

c 2 % 2"' J2 I b k | max | D(e ,u ) _1 | 

* = u 

max | D{e iu )/N{e iu ) | + max | N(e ia y l | 2" £ | a k \ a k 
■-* s »=J (16) 

l-r'ElMA max | N(e { T l I 

4 = a 

2A A Realistic Example 

For the ideally stable second-order bandpass filter 

w n = x n — CLiW n -i — a 2 w n - 2 , n ^ 2 



* It is assumed here that the range of values assigned to the mantissa includes 
the number zero. 
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with poles in the z-plane at angles ^ ±45° and at distance £zl 0.001 
(but not less than 0.001) from the unit circle, we have a l w —1.41, 
a 2 ft* 1, and min u | D(e iu )~ l | » (0.00141)" 1 . We assume that the 
operations are performed as indicated in Fig. 2, so that (3 *» 1, a x — 
3(1.06), and a 2 = 3(1.06). Assuming that t = 27, we find that c our 
bound on the "asymptotic output error-to-signal ratio," ignoring 
input quantization effects, is approximately 0.584 X 10~\ For this 
problem, our bound on c 2 is approximately 0.18 X 10~\ Thus, even 
taking into account input quantization effects, the error-to-signal 
ratio is not more than 0.764 X 10~ 4 . Finally, a simple calculation 
shows that this filter is stable in the presence of roundoff, in the sense 
of inequality (11). 

III. PROOFS 

3.1 Derivation of Inequality (4) 

If a and 6 are floating-point machine numbers, then the floating- 
point product and sum fl(ab) and fl(a + b), respectively, satisfy 

fl(ab) = ab(l 4- e) 

fl(a + b) = (a + 6)(1 + 8) 

with I e I ^ 2" and I 8 I ^ 2"'. Thus, 



\fc-0 fc-l I 



is equal to the value of the output signal 6 of the flow graph described 
in Section II with 

(i) b' k = b k x n _ k 

a' k = dkl/n-k 



b' 




Fig. 2 — Flow graph for the second-order band-pass filter. 
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and 

(it) each of the branch transmissions of the form: (1 + e) with 
| e | ^ 2~ l (recall that in certain special cases e is taken to be 
zero), or -(1 + e) with | e | ^ 2~ ( . Therefore, 

//(£ b&n-u - £ a k y n - k ) 

\i-0 *-l ' 

is equal to 

M N 

J2 btX n - k q h — J2 akt/n-iTk 

A- = fc-1 

in which 

(! _ 2~y w ^ g k ^ (1 + 2~') p ' w (17) 

and 

(I _ 2 -')"-<* J ^ r, ^ (1 + 2- , )" bC * > . (18) 

Inequalities (17) and (18) imply" 

1 - (1.06)^(fc)2 _1 ^ g k £ 1 + (1.06)^)2"' 

1 - (1.06)*> a (/b)2-' ^ r* ^ 1 + (1.0G)^ Q (/,-)2-' 

provided that 2"' max* ^(/c) < 0.1 and 2 - ' max* p a (fc) < 0.1. 
Thus, for n ^ N 

w v 

with 

\vn I ^ 2"' £ I 6, I- 1 *„_* | /3 fc + 2-' Z | a, |- 1 ?/„_, | a k (20) 

and 

fa = (1.06)^(A-), a* = (1.06)p o (A-\ 
Using (1) and (19), 

V 

2] a*e„-* = *?.. , n ^ 



(19) 
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in which, with y„ = w n = for n < 0, 

N 

V, = £ a k (y„- k . - M>„_i) 

for n = 0, 1, • • • , (N — 1). By Propositions 1 and 2 (see Sections 3.5 
and 3.6) 

(e) K ^ max | D(e iu )~ l I (v)k - & ^ 0. (21) 

0S<j<2t 

By Proposition 3 (Section 3.7), inequality (20), and Minkowski's 
inequality 

w^GfrrS 1 " 1 ') 1 

+ 2-' £ I 6 4 | 4t^T7 t I *. I 2 )* +2-'II a. I «*<!/>« (22) 

for all Jv ^ iV. This proves inequality (4). 

3.2 Inequality (8) 

Here we assume that both D(z) and JVC?) are zero free for [ z [ ^ 1, 
that fr ^ 0> and that 

#(«'") | > 2- Z | &* | A . (23) 



m m 

II- u-r'ir t=o 



From (10), we have, with a = 1, 

S a*y.-. = £ &**„-* + 7- - » ^ o- ( 24 ) 

where 

q„ = t/„ , n ^ A r 

= £ atf_* - £ &*_, , n = 0, 1, 2, • • • , (N - 1) 

with .r„ = //„ = for /( < 0. Therefore, by Propositions 1 and 2, 
(x) K ^ max | D(e ia )/N(e ia ) | (i/>k + max | AVT' I (?>« . K ^ 0. 

(25) 
Using Proposition 3, Minkowski's inequality, and (20), 

(q) K ^ (-1-- g | g. I 2 )* + 2"' g | 6. | k (x) K 

+ 2- £) | a k | *&)* , A' ^ iV. (26) 
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Therefore, 

.v 

max | D(e ia )/N(e ia ) | + max | A^e'T 1 | 2"' 2 | o* | a k 



(x)k ^ 



1 - Z" Z I 6* I & max I WT ' | 



<//>. 



,* 



max|JV( e 'T' I (jFXT ^ I ff. I")' 

+ ^ VA + * '-" ^- (27) 



1 - 2~' £ | 6* | ft max | AV")' 
for all 2C ^ AT, which together with (21) and (22) yields 



(e) K ^ 2"' max | Z)(e''")" 



£ I a* | a* + X I b* I 0k 

k—J A— 



max | D(e <M )/N(e im ) | + max | tffc*")' 1 | 2" £ I a* I «* 



<//>* 



1 - 2"' £ | 6* | ft max | AW | 

A = u 

+ max | DC*'")" 1 [ (jc^fj Z I *. 1°) + ™u | Z>(«'T' | 2 



E I &* I ft — 



max I N(e ,a )' 



K + 



T 2 I 2» 

A n = 



1 - 2" 2 1 6 * 1 0* max I #( e ''")" 



(28) 



This proves that there exists a constant c and a function f(K) with 
the property that f(K) — » as K — » oo such that (8) is satisfied for all 
7C ^ A", and of course it also proves that c and f{K) are bounded as 
stated in Section 2.2. 

3.3 Proof of (11) Under the Conditions Stated 
From (24) and Propositions 1 and 2, 

(y) K ^ max | N(e iu )/D(e iu ) | (x) K + max | 0(0" I <?>* ■ 



and using (26) 
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U 

max | N(e ,a )/D(e <u ) | + max | D(e iu y l I 2~* E I h \ & 

(v)k ^ — "— ~ " <*>* 

1 - 2" 1 £ l«* I «* max I D(e-)- | 



^ I DM" I (fVi g I *« 1 2 ) 1 



max 
+ - __ 

1 - 2- E | o 4 | a t max | D^")" 1 | 

which completes the proof. 

3.4 Derivation of Inequalities (15) and (16) 
We have, from (1) and (13), 

E a*(«; B _* - IB..*) = £„ , n ^ (29) 

A (I 

in which a = 1, 

M 

*-o 
and 

£„ - 0, w = 0, 2, ••• , (N - 1). 

Since .f„ = sgn (x n )h2 b for some integer 6 and some /* e [£, 1] (assuming 
that x„ 9^ 0), the magnitude of the error in approximating x n by the 
closest machine number x n = sgn (x n )a2 6 is at most £2~'2 = ^2~'a~ 

| .t„'J ^ 2"' | ;c,. |. Therefore, for n ^ N 

|&|£ 2"« El 6.1-1*-* I. 

1=0 

and by Propositions 1. 2, and 3 
(to - u»)« ^ max | fW | 2" 

■eimGtzt t ix.rY, /c^AT. (so) 

From (30) and (27) 
(w - w) K £ 2~' E | b* | max | D^T | 
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If 

max | D(e iu )/N(e iu ) | + max | A^'T 1 I 2"' £ | a k \ a k 

~ M * -*=*— -{V)K 

1 - 2"' Z I b k | & max | N(e' m ) * | 

2-' £ | ft* | max | DtfT' I max | N(fi im )' 1 I (^TT £ I <?- I') 
1 - 2"' £ I 6, I ft max | A^T' I 

fc=0 u 

for all K ^ A 7 ", provided that N(z) ^ for | z | ^ 1, 6 ** 0, and 
min | N(e iu ) \ > 2~ ( £ | 6* | ft . 

This completes the derivation. 

3.5 Proposition 1: 
If 

I>,r„_, = t^A-i +/«. re ^° 

with: r„ = s„ = for w < 0, c„ ^ 0, and £?-o c t «"' ^ for [ 2 | ^ 1, 
then 

r„ = 2 w »-* s * + £ y »-*/* - n = ° 

23 I W » I < °° I S I y " I < °° . 

tt-0 71 = 

co L' I L 

n-0 1=0 I 1=0 



in which 



and 



for ^ w ^ 2tt. 
Proof:* 

* The proof of this result, although rather trivial, is included because the 
writer knows of no reference where it is proved without the assumption that the 
sequences {s„} and {/„} are 2-transformable. 
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Let M > 0, and let 

§ n = s n for n ^ M 

= for n > M 

J„ = /„ for n ^ M 

= for n > M. 

Then r n = >"„ for n ^ .1/, with 

i. i. ■ 

Z <v\.- 1 = Z rfi«n-i + /-. . n ^ 

and with {r„j, {!„}, and {/„} ^-transformable. Therefore, we have 
in which 

fm = Z v 

,s^) = z Sn2 - 

n = n 

fy) = Z u~ n . 



Thus, 



and hence 



K = Z «.-A + Z "«-*/* , ft ^ 



A- = 11 A=U 



r„ = Z «*.-*** + Z *>»-*/* , (31) 

for ft = 0, 1, • • • , -1/. However, since il/ is arbitrary, (31) is satisfied 
for all ft 2: 0. This proves Proposition 1. 

3.6 Proposition 2: 

It' 

/. = T,c«-i9i . n> 
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with £"-0 | cj | < °° , then 



</>a- ^ max 



Zc 



V 



<<?>, 



for all if ^ 0. 

Proof: 



in which 



Thus, 



„ „ -£7T Jo I n=0 / 'i 

Z7T Jo I n = J-(l 

0, = g t , 1 = 0, 1, ••■ ,K 

= 0, / > /v. 

El/.l'sM I«"" Evrf, 

„ = ^ -'O I n-0 <=<» 

= I- r I ± w-"" t •-*•». 

■S7T Jq | i=0 n-0 



f/co 
rfw 



^ max JZ c ' e 



< max 



E 



c t e 



61T Jq I n-0 

t I 9. I 2 . 



which proves Proposition 2. 
3.7 Proposition 3: 



If 



£ 

I 

1=0 



U\ ^ EU.l-U.-i I, n^iV 



with L ^ N, then 

E 

Vn-AF 

for all /C ^ iV. 



ZI/..I 2 M2U.I S I** 



A' 

c 

n-N-L 
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Proof: 



S I /. I 2 * S 

n-.V n = .V 

K 

& E 



2 I 9i I' I ^--i 

z, I #< r 1 9( II ft--- 



in which 



/,.„ = h n n = 0, 1,2, • •• , /v 

= « > A'. 

Therefore, by the Schwarz inequality, 

t\L\'*tt\9i\ tl^.l-lC 

?i-:V n = .V 1=U 1=0 

/, £ / JC 

^ 2 1 0* I 2 (l £» I 2 I &»-« 

/=U 1=0 \ n=.V 

/. /. / K-l 

^ SI ^1 2 1^1 2 I A. 



(-0 

= (i 1 1, 



I =11 



£ 



This completes the proof. 
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